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Abstract 

It is shown that the nonlocal Dirac operator yielded by a lattice model that preserves 
chiral symmetry and uniqueness of fields, approaches to an ultralocal and invariant under 
translations operator when the size of the lattice tends to zero. 
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The Nielsen-Ninomiya theorem |] has been challenging the construction of quiral field 
theories on the lattice during two decades. The use of Dirac operators D satisfying the 
Ginsparg- Wilson relation || 

75.D + D75 = aDj 5 D, 

is a prospective approach to get chiral fermions with no doublers [chiral symmetry remains 
broken by a term 0(a)] that avoids a contradiction with the Nielsen-Ninomiya theorem 
and yields an exact symmetry of the fermion action ||. However, these operators cannot 
be ultralocal but may be local || (when the interaction range of the operator is bounded 
by some finite lattice distance the operator is termed ultralocal, if it is exponentially 
bounded, the operator is termed local). 

In spite that the locality of the Dirac operator is in general a desired property, nonlocal 
operators have been used to circumvent the restrictions of the Nielsen-Ninomiya theorem 
(see for example || and references there in) and they may have zero modes satisfying the 
Atiyah-Singer index theorem || . This fact strengthen the use of nonlocal operators as an 
approach to follow in lattice QCD. 

Recently, a lattice scheme where the fermion doubling and the chiral symmetry breaking 
are absent has been presented in 0-0. Such a scheme yields a nonlocal Dirac opera- 
tor, but due to some properties of the differentiation matrices used to construct it, the 
quantum algebras and some operational relations maintain their canonical forms on the 
lattice when the number of nodes tends to infinity (or the size of the lattice tends to 
zero) |5|-[rU|. The aim of this letter is to show that such nonlocal differentiation matrices 



and therefore the Dirac operator constructed with them, approach ultralocal operators 
(the nonlocality tends to disappear) when the size of the lattice tends to zero. This fact 
represents another building stone in our lattice model since the chiral symmetry is exactly 
preserved with no doublers and (asymptotically) gauge invariance. 

Differentiation matrices. These matrices, arising naturally in the linear interpolation prob- 
lem of functions, yield exact values for the derivative of polynomial functions at certain 
points selected as lattice sites. Due to this, they have been used to solve some boundary 



value problems (see @-|10[ and references there in). The class of functions on which 
the differential operator to be discretized acts on, determines the kind of differentiation 
matrix to be used. Thus, to get a discrete form of an operator acting on functions that 
fall off rapidly to zero at large distances, we can use the skew-symmetric differentiation 
matrix M 

r o, i = j, 

(D x ) jk = I r_iy+k (1) 
— z — ' i ^ 

constructed with the N zeros Xj of the Hermite polynomial H^(x) as lattice points^, 



whereas to get a discrete form of an operator acting on periodic functions [10], the differ 



1 Equation ([j]) gives a simplified expression for D x . 
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entiation matrix 

', 3 = k, 

(A)ifc = { (-iy+ k . (2) 

2sin^' 

should be used. In this case the lattice points are chosen to be the N equidistant points^ 



tj = -n + 2nj/N, j = 1, 2, ■ ■ ■ , N. (3) 



Since the fields are expected to fall off rapidly to zero at large distances or one may 
assume periodic boundary conditions on the fields imposed by the lattice, we can use 
differentiation matrices of the type ([I]) or ([J), or a combination of them, to construct 
a discretized formulation of classical field theory in 3 + 1 variables where the fields are 
unique (no doublers), the chiral symmetry is not broken when m = and the gauge 
invariance is recovered in the limit iV — > oo. Some of these properties have been shown 
previously 0-f| for the differentiation matrix (jH). 

In this letter we analyze the ultralocality property of these differentiation schemes used to 
construct the Dirac operator. In addition, we show the absence of the fermion doubling 
in the formulation with derivative of the type (|2|). 

First of all, let us show the relationship between the differentiation matrices D x and D t . 
Let j and A; be a pair of distinct and fixed integers. Since the nth zero of the Hermite 
polynomial H N (x) can be approximated by 

s.- ' 2 ":^ 1 '* <4) 



2V2N 

when n is fixed and N is large enough, we have that 

/ 2N(-lV +k 



D x ) jk = 2J2/N(D t ) jk = (-iy +k 



Ti j — k 



This relation does not become a matrix equation between D x and D t since is valid 
only for fixed n when iV — > oo and therefore, sin(tj —tk)/2 can not be approximated by 
(tj — tk)/2 for elements close to the border of D t . Thus, if iV is large enough, 



D x = 2p/ND t + O, (5) 

where O is a border matrix whose nonzero elements depend on N. By a N x N border 
matrix we mean a matrix M whose elements Mj k lying between n < j, k < N — n are all 
of them zero for certain n « N . 



Clearly, D t is a Toeplitz matrix [] L] and fulfills the requirements of translational invariance 



on the periodic lattice (^|), whereas D x becomes a Toeplitz matrix (except for a border 
2 The same result is obtained for the set of points tj = jf(2j — N — 1). 
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matrix) only for iV large enough. Both of them are nonlocal matrices for a finite N 
since | csc(rc) | > l/\x\ > exp(— \x\) in the corresponding ranges. However, we can use 
the Toeplitz property to show that both D x and D t approach ultralocal matrices almost 
everywhere when N — > oo. By almost everywhere we mean except for border matrices. 
Such border matrices circumvent the no-go theorem given in [|J and allow us to have an 
ultralocal Dirac operator almost everywhere with the properties given above. To show 
this, we will use the matrix D t which has a stronger nonlocality. 



Ultralocality of D t . Since matrix D t is a skew-symmetric Toeplitz matrix, it can be 
written as 

Af-l 

D t = Y,d n (A n -B n ), (6) 

n=l 

where 

(_l)n+i 

dn = 2sin(nvr/iV) (7) 
and A and B are the backward shift and forward shift matrices respectively, i.e, 

Ajk = Sj t k-i, Bjk = 

These matrices satisfy the useful relations 

ABA = A, BAB = B, (8) 

which can be used to show that 

»n— 1 m— 1 



A 2m _ B 2m = _ 5 ) ^ (^2n+l + B 2n + 1^ _ ^ B](l + ^ + B 



_j2n 

n=0 n=l 

and 

m rn—l 



A 2m+i_ B 2 m +i = (A-B)(l+J2(A 2n + B 2n )) - [A,B] J2(A 2n+1 + B 



j2n+l _,_ r2ti+1\ 

n=l n=0 



Here, I stands for the N x N identity matrix. Thus, @ can be rewritten in the useful 
form 

D t = (A-B)A-[A,B]E, (9) 
where A and £ are the symmetric matrices 

N-2 N-3 

A = s 1 I+Y / s n+1 (A n + B n ), Z = s 2 I+Y, s n+2{A n + B n ), (10) 

n=l n=l 

and 

[N/2]-l [(AT-l)/2] 

E d 2n . (11) 



Here, [x] stands for the greatest integer less or equal to x. Additionally to (§), A and B 
satisfy 

1, j = k = l, 



([AB]) 



J 



k = N, 



0, otherwise. 
Therefore, the product 

1 = [A,B]Z 

appearing in (§) becomes a border matrix. For N great enough, the sums (|ll]) appearing 
in (|1(]) can be approximated by 



Ne n ,,n7i\ 

^T cot( 2iV' 



n 



l,2,...,iV-l, 



where 



1, N + n even, 

2, N + n odd. 



Since D t = (A — i?)A (except for the border matrix Ox), the interaction range of D t is 
essentially the one of A, which is measured by |su_fc + i||. If we put x = rrn/N, \s n \ can be 
rewritten as 

Ae "cot(x/2) 



six 



N 

47T 



log 



7T 



x G (0,7r). 



Note that s(x) is an antisymmetric function in (0,7r) with respect to 7r/2 and that the 
absolute values of Si = s(0) and s^-i = s(7r) grow as (N/Ati) \og(2e n N 2 / it 2 ) . Thus, for x 
other than or ir, 

S{X) - - ~ (12) 



s(0) 



0. iV^oo. 



This is the behavior of an unnormalized Dirac's delta. Thus, s(x) — > Cn[S{x) + 5(x — 7r)]. 
To get the constant of normalization the measure <5x = 7r/iV has to be taken into account. 
For convenience, we integrate from —tt/2 to tt/2 to obtain 



N\ ,Ne r 
— log( 

47T 7T 



-1 



tt/2 



s(x)dx 



1. 



(13) 



Eqs. ([12] ) and (|I~3|) show that the asymptotic behavior of s(x), a; G (— tt/2, 7r/2), is that 
of the Dirac delta multiplied by a growing numerical factor. Therefore, except for the 
measure, 



\A(x-y)\ -> ^log(— ^)[5(a; - y) +5(|s-?/| 

47T 7T L 



7T 



N 



oo. 



(14) 



Since s(7r) corresponds to the elements of A located on the band defined by sn-i, the 
second delta of the right-hand of ([14]) yields another border matrix. Let us denote by 
O2 such a matrix. Then, we have from @ that D t approaches an ultralocal matrix 
as iV — > 00 except for the sum of the border matrices 0\ and O2. If D t is applied to a 
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vector formed with the values of a function f(x), the matrix- vector multiplication changes 
into an integration in the asymptotic limit whenever the measure is taken into account. 
Therefore, the action of D t on such a vector is mainly given by 

£ D t (x - y)f(y) - [f(x + a) - /(x - o)]/2o, a ^ 0, (15) 

where 

-i Ne n 

a =4^ log( ^T ) - 

This is an expected result since D t is an exact differentiation matrix for certain subclasses 
of periodic functions UlO |. Due to @, the differentiation matrix D x also satisfies a relation 
like QT3). 



Uniqueness of the free fermionic fields. Next, we show that the use of the differentiation 
matrix (0) avoids the doubling fermion problem on a lattice. Let us consider a four- 
dimensional regular lattice with N 4 points whose support are in the interval — L < x% < L 
[substitute tt by L in ([$])]. We assume functional periodicity in the fermion fields as 
trigonometric polynomials. The (continuum) Euclidean free Lagrangian reads 

£ = ip(x) (7^ + m) ip(x) . 

We make the substitution 

W-n M n M — 1 

where we have used (El) and (|7|), and /x = 1, 2, 3, 4. Thus the discretized action is 



(16) 



s = J2 



n,fi 



N-n„ 



n u — 1 



lu=l 



la=l 



At first glance this is a rather nonlocal action, since the field at each lattice point couples 
with the field evaluated at the reminder points along each axis. Nevertheless, due to the 
fact that [D^Jn is the exact partial derivative along the direction ji at the point n (see 



field 



IQ[), the action can be easily diagonalized by taking the discrete Fourier transform in the 

1 

7n* , 



and in its derivative 



(17) 



here = nl^/L, with Z M being an integer restricted to values in the range — (N — l)/2 < 
I < {N — l)/2 (N is assumed to be an odd number). Thus the maximum momentum 
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allowed is \k max \ = n/2A, with A the length between two successive points. It is clear 
that the values of k become to be with no restriction in the limit when N goes to infinity 
and sum in (|17]) must be replaced by an integral in the corresponding Brillouin zone. In 
fact, in this limit the derivative in Eq. ( ]TB| ) tends to the exact continuum value, and the 
continuum theory is recovered. The Fourier transform of the action can be written as 

which yields to the dispersion relation 

S{k) = 1 ^ + m. (18) 

It is clear that no doublers appear since the values of k are proportional to an integer. 
Moreover in the limit in which the spacing A tends to zero (|T8|) gives the exact energy 
values. 

Other desired properties, as gauge invariance and chiral symmetry, also have support on 
this scheme. They can be shown as in Ref. ||. Let us emphasize that the exact chiral 
symmetry is preserved explicitly in our scheme as it happens with the so-called SLAC 
derivative | 12| , which is defined through an equation like fll7|). In that context, matrix 



is an explicit realization of such a derivative with a number of properties given in JTU 
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